Local asymptotic normality for qubit states 
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We consider n identically prepared qubits and study the asymptotic properties of the joint state 
p®n yj e s j low that for all individual states p situated in a local neighborhood of size 1 / y/n of a fixed 
state p , the joint state converges to a displaced thermal equilibrium state of a quantum harmònic 
oscillator. The precise meaning of the convergence is that there exist physical transformations 
T n (trace preserving quantum channels) which map the qubits states asymptotically close to their 
corresponding oscillator state, uniformly over all states in the local neighborhood. 

A few consequences of the main result are derived. We show that the optimal joint measurement 
in the Bayesian set-up is also optimal within the pointwise approach. Moreover, this measurement 
converges to the heterodyne measurement which is the optimal joint measurement of position and 
momentum for the quantum oscillator. A problem of local state discrimination is solved using local 
asymptotic normality. 

PACS numbers: 



1. INTRODUCTION 

Quantum measurement theory brings together the 
quantum world of wave functions and incompatible ob- 
servables with the classical world of random phenomena 
studied in probability and statistics. These fields have 
come ever closer due to the technological advances mak- 
ing it possible to perform measurements on individual 
quantum systems. Indeed, the engineering of a novel 
quantum state is typically accompanied by a verification 
procedure through which the state, or some aspect of it, 
is reconstructed from measurement data 0. 

An important example of such a technique is that of 
quantum homodyne tomography in quantum òptics 0- 
This allows the estimation with arbitrary precision of the 
whole density matrix [Ü Ql ü, El of a monochromatic 
beam of light by repeatedly measuring a sufficiently large 
number of identically prepared beams Q, 0, • 

In contrast to this "semi-classical" situation in which 
one fixed measurement is performed repeatedly on inde- 
pendent systems, the state estimation problem becomes 
more "quantum" if one is allowed to consider joint mea- 
surements on n identically prepared systems with joint 
state p® n . It is known j^] that in the case of unknown 
mixed states p, joint measurements perform strictly bet- 
ter than separate measurements in the sense that the 
asymptotical convergence rate of the optimal estimator 
p n to p goes in both case as C ' / 1 \fn with a strictly smaller 
constant C in the case of joint measurements. 

Let us look at this problem in more detail: we dispose 
of a number of n copies of an unknown state p and the 
task is to estimate p as well as possible. The first step 
is to specify a cost function d(p n , p) which quantifies the 
deviation of the estimator p n from the true state. Then 
one tries to devise a measurement and an estimator which 
minimizes the mean cost or risk in statistics jargon: 

R{p,Pn) ■■= (d(p n (X),p)), 



with the average taken over the measurement results X. 
Since this quantity still depends on the unknown state 
one may choose a Bayesian approach and try to optimizc 
the average risk with respect to some prior distribution 
7T over the states 
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R(p,p n )ir(dp). 



Results of this type have been obtained in both the pure 
state case E(l Ull E2 E3 El EM EJ an d the mixed 
state case [Ï8t IÏ9t l20l Í2lL l22l l23l l24j . However most of 
these papers use methods of group theory which depend 
in on the symmetry of the prior distribution and the form 
of the cost function, and thus cannot be extended to ar- 
bitrary priors. 

In the pointwise approch [ÏJ Hj||2(||23, 153, |2íJ one tries 
to minimize R(p, p n ) for each fixed p. We can argue that 
even for a completely unknown state, as n becomes large 
the problem ceases to be global and becomes a local one 
as the error in estimating the state parameters is of the 
order ^= . For this reason it makes sense to parametrize 
the state as p := p{&) with belonging to some set in 
R fc and to replace the original cost with its quadratic 
approximation at 0: 

d(0,è n ) = (9-è n ) T G(9)(8-Ò n ), 

where G is a kx k positive, real symmetric weight matrix. 

Although sccmingly different, the two approaches 
can be compared 30], and in fact for large n the 
prior distribution ir of the Bayesian approach should 
become increasingly irrelevant and the optimal Bayesian 
estimator should be close to the maximum likelihood 
estimator. An instance of this asymptotic equivalence is 
proven in Subsection l7 Al 

In this paper we change the perspective and instead 
of trying to devise optimal measurements and estimators 
for a particular statistical problem, we concentrate our 



2 



attention on the family of joint states p(8)® n which is 
the primary "carner" of statistical information about 6. 
As suggested by the locality argument sketched above, 
we consider a neighborhood of size -4= around a fixed 
but arbitrary parametcr 9q, whose points can be written 
as 8 = 6*o + u/y/n with u £ R fc the "local parameter" 
obtained by zooming into the smaller and smaller balls 
by a factor of \fn. Very shortly, the principle of local 
asymptotic normality says that for large n the local family 



Pn ■= P (0o + u/Vn) 



<C, 



converges to a family of displaced Gaussian states <fi u of 
a of a quantum system consisting of a number of couplcd 
quantum and classical harmònic oscillators. 

The term local asymptotic normality comes from 
mathematical statistics j3Ï| where the following re- 
sult holds. We are given independent variables 
X\ , . . . , X n £ X drawn from the same probability dis- 
tribution p e «+ u / vn over % depending smoothly on the 
unknown parameter u £ R fc . Then the statistical infor- 
mation contained in our data is asymptotically identical 
with the information contained in a single normally dis- 
tributed Y £ R fe with mean u and variance I(#o) -1 > the 
inverse Fisher information matrix. This means that for 
any statistical problem we can replace the original data 
... , X n £ X by the simpler Gaussian one Y with the 
same asymptotic results! 

For the sake of clarity let us consider the case of 
qubits with states parametrized by their Bloch vectors 
p(r) = |(1 + r cr) where <x = {o~ x , a y , a z ) are the 
Pauli matrices. Define now the two-dimensional family 
of identical spin states obtained by rotating the Bloch 
vector ro = (0, 0, 2/z — 1) around an axis in the x-y plane 



U 



H 
1-// 



u 



u £ 



(l·l) 

with unitary Í7(v) := cxp(i(v x cr x +v y o~ y )) and „ < fi < 1. 

Consider now a quantum harmònic oscillator with po- 
sition and momentum operators Q and P on L 2 (R) satis- 
fying the commutation relations [Q,P] = il. We denote 
by {\n),n > 0} the eigenbasis of the number operator 
and define the thermal equilibrium state 



0°=(i-i>)xyiA><*i, 

fc=0 

where p = - zií . We translate the state 0° by using the 
displacement operators D{z) = exp(za* — za) with z £ C 
which map the ground state |0) into the coherent state 



4> u := D(y/2p~ïa u )0 o D{y/2p~ïa u )* , (1.2) 
where a u := — u y + iu x . 



Theorem 1.1 Let p" 6e the family of states i|l.lfl on the 

Hilbert space (C 2 ) ™ and ip^ the family (|1.2|l of displaced 
thermal equilibrium states of a quantum oscillator. Then 
for each n there exist quantum channels (trace preserving 
CP maps) 



(1.3) 



T n :M[(cy n )^T(L 2 (R)) 
with T(L 2 (M)) the trace-class operators, such that 



hm sup H^-Ï^Oí») !!! = (), 

lim sup ||p^-5„(0 u )|| 1 = O. 
n ^°° ue/ 2 



(1.4) 



for an arbitrary bounded interval I C R. 

Let us make a few comments on the significance of the 
above result. 

i) The "convergence" JEU of the qubit states holds in 
a strong way (uniformly in u) with direct statistical and 
physical interpretation. Indeed the channels T n and S n 
represent physical transformations which are anàlogues 
of randomizations of classical data [3Ï| . The meaning of 
(|1.4|) is that the two quantum models are asymptotically 
equivalent from a statistical point of view. 

ii) Indeed for any measurement M on i 2 (R) we can con- 
struct the measurement M oT n on the spin states by first 
mapping them to the oscillator space and then perform- 
ing M. Then the optimal solution of any statistical prob- 
lem concerning the states p" can be obtained by solving 
the same problem for (f> u and pulling back the optimal 
measurement M as above. We illustrate this in Section 
0for the estimation problem and for hypothesis testing. 

iii) The proposed technique may be useful for appli- 
cations in the domain of coherent spin states |32( and 
squeezed spin states j33|. Indeed, it has been known 
since Dyson [34j that n spin-1 partides prepared in 
the spin up state | I)* 8 ™ behave asymptotically as the 
ground state of a quantum oscillator when considering 
the fluctuations of properly normalized total spin com- 
ponents in the directions orthogonal to z. Our Theorem 
extends this to spin directions making an "angle" u/y/n 
with the z axis, as well as to mixed states, and gives a 
quantitative expression to heurístic pictures common in 
the physics literature (see Section |ÏÏJ). We believe that 
a similar approach can be followed in the case of spin 
squeezed states and continuous time measurements with 
feedback control [35| . 

Next Section gives an introduction to the statistical 
ideas motivating our work. In Section[3]we give a heurís- 
tic picture of our main result based on a the total spin 
vector representation of spin coherent states familiar in 
the physics literature. 
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The proof of Theorem 11.11 extends over the Sections 
I4I5I6I and uses methods of group theory and some ideas 
from [Ü HE HI 

Section describes a few applications of our main re- 
sult. In Subsection l7 Al we compute the local asymptotic 
minimax risk for the statistical problem of qubit state 
cstimation. An estimation scheme which achieves this 
risk asymptotically is optimal in the pointwise approach. 
We show that this figure of mèrit coincides with the risk 
of the heterodyne measurement and that it is achieved 
by the optimal Bayesian measurement for the SU (2)- 
invariant prior pü . l2flj . This proves the asymptotic 
equivalence of the Bayesian and pointwise approaches. 

In Subsection 17 Bl we continue the investigation of the 
optimal Bayesian measurement and show that it con- 
verges locally to the heterodyne measurement on the os- 
cillator which is a op timal joint measurement of position 
and momentum [39| . 

Another application is the problem discriminating be- 
tween two states p„ which asymptotically converge to 
each other at rate \j\fn. In this case the optimal mea- 
surement for the parameter u is not optimal for the test- 
ing problem, showing in particular that the quantum 
Fisher information in general does not encode all sta- 
tistical information. 



2. LOCAL ASYMPTOTIC NORMALITY IN 
STATISTICS AND ITS EXTENSION TO 
QUANTUM MECHANICS 

In this Section we introduce some statistical ideas 
which provide the motivation for deriving the main re- 
sult. 

Quantum statistical problems can be seen as a game 
between a statistician or physicist in our case, and Na- 
ture. The latter tries to codify some information by 
preparing a quantum system in a state which depends on 
some parameter u unknown to the former. The physicist 
tries to guess the value of the parameter by devising mea- 
surements and estimators which work well for all choices 
of parameters that Nature may make. In a Bayesian set- 
up Nature may build her strategy by randomly choosing 
a state with some prior distribution. In order to solve the 
problem the physicist is allowed to use the laws of quan- 
tum physics as well as those of classical stochastics and 
statistical inference. In particular he may transform the 
quantum state by applying an arbitrary quantum channel 
T and obtain a new family T(p u ). In general such trans- 
formation goes with a loss of information so one should 
have a good reason to do it but there are non trivial sit- 
uations when no such loss oceurs |40| , that is when there 
exists a channel S which reverses the effect of T restricted 
to the states of interest S(T(p u )) = p u . If this is the case 
the we consider the two families of states p u and T(p u ) 
as statistically equivalent. 

In statistics such transformations are called random- 
izations and a useful particular example is a statistic. 



which is just a function of the data which we want to an- 
alyze. When this statistic contains all information about 
the unknown parameter we say that it is sumeient, be- 
cause knowing the value of this statistic alone suffices and 
given this information, the rest of the data is useless. For 
example if Xi, . . . X n £ {0, 1} are results of independent 
coin tosses with a biased coin, then X — — V . Xi is suffi- 
cient statistic and may be used for any statistical decision 
without loss of efficiency. 

Quantum randomizations through quantum channcls 
allows us to compare sccmingly different families of states 
and thus opens the possibility of solving a particular 
problem by casting it in a more familiar setting. The 
example of this paper is that of state estimation for n 
identical copies of a state which can be cast asymptoti- 
cally into the problem of estimating the center o a quan- 
tum Gaussian which has a rather simple solution |39| . 
The term "asymptotically" means that for large n we 
can find quantum channels T ni S n which almost map the 
the families of states into each other as in equation (|1.4fl . 

The second main idea that we want to introduce is 
that of local asymptotic normality. Back in the coin 
toss example we have that X is a good estimator of the 
probability p of obtaining a 1 and by the Central Limit 
Theorem the error X — p has asymptotically a Gaussian 
distribution 

VE(X-p)^N(0,l/p(l-p)), 

in particular the mean error is ((X— p) 2 ) = l/(np(l—p)). 
Now, if for each n the unknown parameter p is restricted 
to a local neighborhood of a fixed po of size l/\/n, one 
might expect an improvement in the error because we 
know more about the parameter and we can use that 
information to built better estimators. However this is 
not entirely true. Indeed if we write p = po + u/ -Jn then 
the estimator of the local parameter u is 

ü n = \fn{X - po) N(u, l//io(l - po)) 

which says that the problem of estimating p in the lo- 
cal parameter model is as dimcult as the original prob- 
lem, i.e. the variance of the estimator is the same. The 
reason for this is that the additional information about 
the location of the parameter is nothing new as we could 
guess that dircctly form the data with very high probabil- 
ity. Thus without changing the difhculty of the original 
problem we can look at it locally and then we see that it 
transforms into that of estimating the center of a Gaus- 
sian with fixed variance N(u,l/ po(l — po)), which is a 
classical statistical problem. 

In general we can formulate the following principle: 
given X\ , . . . , X n £ X independent with distribution 
p0 o +u/Vn depending smoothly on the unknown parame- 
ter u £ K fe , then asymptotically this model is statistically 
equivalent (there exist explicit randomizations in both 
directions) with that of a single draw Y £ R k from the 
Gaussian distribution iV(u, I(6 )~ 1 ) with fixed variance 
equal to the inverse of the Fisher information matrix |3l| . 



In the quantum case we replace the randomizations 
by quantum channels and the Gaussian limit model by 
its quantum equivalent which in the simplest case is a 
family of displaced thermal states of a quantum oscillator 
(see Theorem ll.lfl . but in general is a Gaussian state on 
a number of coupled quantum and classical oscillators, 
with canonical variables satisfying general commutation 
relations |4l| . 

A simple extension of Theorem 11.11 is obtained by 
adding an additional local parameter t <E R for the den- 
sity matrix eigenvalues such that fi — fiQ + tj y/rï. This 
leads to a Gaussian limit model in which we are given a 
quantum oscillator is in state 4> u and additionally, a clas- 
sical Gaussian variable with distribution N(t, l//io(l — 
/io))- The meaning of this quantum-classical coupling 
is the following: asymptotically the problem of estimat- 
ing the eigenvalues decouples from that of estimating the 
direction of the Bloch vector and becomes a classical sta- 
tistical problem (identical with the coin toss discussed 
above), while that of estimating the direction remains 
quantum and converges to the estimation of a Gaussian 
state of a quantum oscillator. We note that this decou- 
pling has been also observed in 0, HÜ ■ 



3. THE BIG BALL PICTURE OF COHERENT 
SPIN STATES 

In this section we give a heuristic argument for why 
Theorem II . II holds which will guide our intuition in later 
computations. 

It is customary to represent the state of two dimen- 
sional quantum system by a vector ~~r in the Bloch sphere 
such that the corresponding density matrix is 



1, 

P= r a > 



1 



(1 + r x a x + r y (7 y + r z a z ), 



where <Ji represent the Pauli matrices and satisfy the 
commutation relations [<Ji,<jj] — lie^Ok. In particu- 
lar if r = (0,0, ±1) then the state is given by the spin 
up 1 1) and respectively spin down | J.) basis vectors of C 2 , 
and the z-component of the spin <j z takes value ±1. As 
for the x and y spin components, each one may take the 
vàlues ±1 with equal probabilities such that on average 
(ox) = (&y) = but the variances are (al) = (a y ) = 1. 
Moreover a x and a y do not commute and thus cannot be 
measured simultaneously. 

What happens with the Bloch sphere picture when we 
have more spins? Consider for the beginning n identical 
spins prepared in a coherent spin up state | T)® n j then we 
can think of the whole as a singlc spin system and define 

the global observables L^' = J2k=i a i f° r * *= x ^V^ z y 

(k) 

where er) is the spin component in the direction i of 
the fc's spin. Intuitivcly, we can represent the joint state 
by a vector of length n pointing to the north pole of a 
large sphere as in Figure^ However due to the quantum 
character of the spin observables, the x and y components 




FIG. 1: (Color online) Quasiclassical representation of n spin 
up qubits 



cannot be equal to zero and it is more instructive to think 
in terms of a vector whose tip lies on a small blob of the 
size of the uncertainties in x and y, sitting on the top of 
the sphere. Exactly how large is this blob ? By using 
the Central Limit Theorem we conclude that in the limit 
n — > oo the distribution of the "fiuctuation operator" 



c(ra) 



1 



>2n 



1 



converges to a N(0, 1/2) Gaussian, that is (S x ) = and 
(S x ) w 1/2, and similarly for the component S y l \ The 
width of the blob is thus of the order y/n in both x and 
y directions. 

Now, the two fluctuations do not commute with each 
other 



(3T) 



which is the well know commutation relation for canoni- 
cal variables of the quantum oscillator. In fact the quan- 
tum extension of the Central Limit Theorem |36l | makes 
this more precise 



lim 

n — >oo 



<8m / 



(t i n < i tr i = (fi. n x ^ fi)- v ^ e f>. 



k=l 



k=l 



where X x := Q and X y := P satisfy [Q, P] = il and 
is the ground state of the oscillator. 

The above description is not new in physics and goe s 
back to Dyson's theory of spin-wave interaction |34| . 
More recently squeezed spin states |23] for which the vari- 
ances (S x ) and (S y ) of spin variables are different have 
been found to have important applications various fields 
such as magnetometry |35j . entanglement between many 
partides |42| . The connection with such applications will 
be discussed in more detail in Section 

We now rotate all spins by the same small angle for 
each particle as in Figure [21 As we will see, it makes 
sense to scale the angle by the factor -7= i.e. to consider 



CXp ( -y=(u x a x +U y CTy) ) It) 



u e 
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is a thermal equilibrium state 

oo 

<t>° = (l-p)J2p k \k)(k\, 



P = 



Finally the rotation by exp í-^=(u x a x + u y a y )j produces 
a displaccmcnt of the thermal state such that (Q) = 
-V2(2/x - í)u y and (P) = \/2(2/i - l)u x . 



FIG. 2: (color online)Rotated coherent state of n qubits 



Indeed for such angles the z component of the vector will 
change by a small quantity of the order y/n < n so the 
commutation relations (|H . ip remain the same, while the 
uncertainty blob will just shift its center such that the 
new averages of the renormalized spin components are 
(Sx) ps — y/2u y and (Sy) ~ V%u x . All in all, the spins 
state converges to the coherent state \a u ) of the oscillator 
where ev u = (— u y + iu x ) G C and in general 



\a) :=e W {-\a\ 2 /2)Y,^\3), 

3=0 V3 ' 

with \j) representing the j's energy level. 

We consider now the case of qubits in individual mixed 
state fjt\ TXT | + DU I with < V 2 M < !• Thenthe 

"length" of L z is n(2p, — 1) but the size of the blob is the 
same (see Figurc[3J). However the commutation relations 



(2/1- l)n 



4. LOCAL ASYMPTOTIC NORMALITY FOR 
MIXED QUBIT STATES 

We give now a rigorous formulation of the heurístics 
presented in the previous Section. Lct 



u 

1-fJL 



(4.1) 



be a density matrix on C 2 with /i > 1/2, represent- 
ing a mixture of spin up and spin down states, and 
for every u = (u x ,u y ) £ M 2 consider the state p u = 
C/(u)p°C/(u)*where 



U(u) := exp(i(u x cr x +Uy<7 y )) 



sm u 



-up , 



cos u 



with ip — Axg(—u y + iu x ). We are interested in the 
asymptotic behavior as n — > oo of the family 



, U G J 



(4.2) 



where / = [— a, a] is a fixed finitc interval. 

The main result is that T n is asymptotically normal, 
meaning that it converges as n — > oo to a limit family 
G n '■= {0 u ,u G I 2 } of Gaussian states of a quantum 
oscillator with creation and annihilation operators satis- 
fying [a, a*] = 1. Let 



:= (l-p)J2p k \k){k\, 



(4.3) 



k=0 



be a thermal equilibrium state with \k) denoting the k's 
energy level of the oscillator and p — ——^ < 1. For every 
u G I 2 definc 



FIG. 3: (color online)Quasiclassical representation of n qubit 
mixed states 



of S x and S y do not reproduce those of the harmònic 
oscillator and we need to renormalize the spin as 



5' 



(n) 



1 



v/2(2/i- l)n 



1 



^2(2^ - ïjn 



--L,. 



The limit state will be a Gaussian state of the quantum 
oscillator with variance (Q 2 ) = (P 2 ) = 2 (2^-i) < 3' 



:= D(^2u - la^M-^fi - ïc 



(4.4) 



where D(z) := exp(za* — z*a) is the displacement op- 
erator, mapping the vacuum vector |0) to the coherent 
vector |z) and a u = {—u y + iu x ) . 

The exact formulation of the convergence is given in 
Theorem 11.11 Thus the state p" of the n qubits which 
depends on the unknown parameter u can be manipu- 
lated by applying a quantum channel T n such that its 
image converges to the Gaussian state 4> u , uniformly in 
u G I 2 . Conversely by using the channel S n , the state 
can be mapped to a joint state of n qubits which is 
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converges to p" uniformly in u G I 2 . By Stinespring's 
thcorem we know that the channels are of the form 

T(p) = Trjç (VpV*), 
S(</>) = Trv (yV4>W*) , 

where the partial traces are taken over some ancillary 
Hilbert spaces /C,/C' and 



V : (C 2 f " 
W : L 2 (R) -» 



L 2 (M) ®/C, 



are isometries (V*V — 1 and W*W — 1). 

Our task is now to identify the isometries V n and W„ 
implementing the channels T n and respectively S n satis- 
fying 1)1.4(1 . The first step towards identifying these V n is 
to use group representations methods so as to partially 
(block) diagonalize all the p" simultaneously. 



A. Block decomposition 

In this Subsection we show that the states p" have a 
block-diagonal form given by the decomposition of the 

space (C 2 )®" into irreducible representations of the rel- 
evant symmetry groups. The main point is that for large 
n the weights of the different blocks concentrate around 
the representation with total spin j„ = n(fi — 1/2) . 

The space (C 2 )**™ carries a unitary representation n n 
of the one spin symmetry group SU(2) with n n (u) = 
u ®n f or an y u g 5(7(2), and a unitary representation of 
the symmetric group S(n) given by the permutation of 
factors 



7r„(r) : «i®- • -®w„ h-> u t - i(i)( 



•'T-·Mra) ) 



r G S(n). 



As [n n (u), 7r„(r)] = for all u G SU(2),t G 5(n) we have 
the decomposition 



n/2 
5=0,1/2 



(4.5) 



where the direct sum runs over all positive (half )-integers 
j up to n/2, and for each fixed j, Hj — C 2j+1 is a 
irreducible representation of SU(2) with total angular 
momentum J 2 = j(j + 1), and W n = C rij is the irre- 
ducible representation of the symmetric group S(n) with 
rij = ( n /2_j) ~ \ n /2—j-i) ■ I n particular the density ma- 
trix p" is invariant under permutations and can be de- 
composed as a mixture of "block" density matrices 



n/2 

Pn= ® Pn(j)pln®—, 
5=0,1/2 7 

with probability distribution p n (j) given by j24[: 



(4.6) 



Pn(j) 



2fl - 



T (l-ri T " J (4.7) 



where p := — rp- A key observation is that for large n and 
in the relevant range of j's, p n (j) is essentially a binomial 
distribution 



o,. 



Indeed we can rewrite p n (j) as 

P«Ü) := B nil ,(n/2+j) x K(j,n,n) (4.8) 
where the factor ií(j, n, p) is given by 

n + (2Ü-Ín) + l)/(2/í-l) 



K(j,n,p) := (1-p 2j+1 ) 



«■ + (Í ~jn + 1)/M 



and j n := n(/i — 1/2). As -B n ,u is the distribution of 
the sum of n independent Bernoulli variables with indi- 
vidual distribution (1 — p, p) over {0, 1}, we can use the 
central limit Theorem to conclude that its mass concen- 
trates around the average fin with a width of order i/n, 
in other words of any < e < 1/2 we have 



,1/2+. 



lim 

n — >oo 



B n,u(fJ.n + p) = 1. 



(4.9) 



Let us denote by J n ^ the set of vàlues j of the total 
angular momentum of n qubits which lic in the interval 
[in — rt 1 / 2+c , j„ + n 1 / 2+c ]. Thcn for large rc, the factor 
K(j, n, fi) is close to 1 uniformly over j G J n ,t and from 
formulas 1(4.8(1 . 1(4.9(1 we conclude that p n (j) concentrates 
asymptotically in an interval of order n 1 ^ 2+e around j n : 



lim p n (J n ,e) = 1. 



(4.10) 



This justihes the big ball picture used in the previous 
section. 



B. Irreducible representations of SU(2) 

Here we remind the reader some details about the irre- 
ducible representation nj of SU(2) on Hj. Let a x , a y , a z 
be the Pauli matrices and denote Jjj := TTj(cri) for 
l = x,y,z the generators of rotations in the irreducible 
representation 7ij , such that the corresponding unitàries 
are Uj(u) :— exp (i(u x Jj^ x + u y Jj^y)). There exists an or- 
thonormal basis m), m — —j, . . . , j} of Tíj such that 

Jj,z\j' m ) = m \j, m )- 

Moreover, with Jj_± :— Jj X ± iJj, y we have 



Jj,+ \í- m ) = - m \ / j + m + l \j, m + 1), 

J j-\h m ) = Vj - m + 1 Vj + m li) m - !)• 

With these notations and p = as before, the state 
p° can be written as 29] 

3 

m=—j 
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where the normalizing factor is Cj (p) = (1— p)/(l— p 2j+1 ). 
The rotated block states can be obtained by applying the 
unitary transformation 

Pi n = u j (u/v^)pi n u j (u/v^r, 

with Uj(u) as above. Finally, we define the vectors 



|j,w) := Uj(w)\j,j) 



(4.11) 



which will be used in later computations, and notice that 
their coordinates with respect to the \ j, m) basis are given 
by 0: 



0-'"l7·w) = 1 /( - m jC'-'"(1 IC| 2 )^ (4.12-1 



2 j Vi-mn 
+ m) 

where £ = e liPw sin |w| with ^ = Arg(— Wy + iw x ). 
5. CONSTRUCTION OF THE CHANNELS T n 



Proposition 5.1 The following uniform convergence 
holds 

Hm sup max ||V^ n V;* - u ||i = 0. 

where J n ,t is the set defined above equation (|4.1U|I . 

The proof of the Proposition requires a few ingredients 
which in our opinion are important on their own for which 
reason we formulate them apart and refer to relevant pa- 
pers for the proofs. 

Theorem 5.2 '36] Let a,b E M(C d ), satisfying Tr(a) = 
Tr(6) = and define 

L(a, b) — exp(io) exp(iò) — exp(ia + ib) exp ( — [a, b] j . 



Oi 



we define the fluctuation operator 
F n (a) = -^E a *' 



For each irreducible representation space Tij we define 
the isometry Vj : Hj -> i 2 (M) by 



Vj : |i, m) i-> | j - m) 



(5.1) 



where {|n),n>0} represents the energy eigenbasis of 
the quantum oscillator with eigenfunctions ip n (x) = 



/n2 n n\ E L 2 {R). Using the decompo- 
sition 1|4.5[) we put together the different blocks we con- 
struct for each n E N the "global" isometry 

n/2 n/2 

V n := Vj®l: Hj <8> C n * L 2 (R) ® £ n , 

3=0,1/2 3=0,1/2 

where /C„ := 0™f.Q x / 2 C™?. By tracing over fC n we ob- 
tain the channel T n (p) := Tric n (V n pV£) mapping a joint 
state of n spins into a state of the quantum oscillator. 
This channel satisfies the convergence condition l|1.4fl as 
shown by the estimate 



n/2 



n/2 

E PnÜWjPljV? 
3=0,1/2 



< E p»U)\\VjpïjV? - u IL 

3=0,1/2 

<2 E P«0')+sup max Wjp^V* - ^h, 

where the first term on the right side converges to by 
Q4.10[) . and for the second one we apply the following 
ProDOSition l5.1l which is the major technical contrubition 
of this paper. 



where m = l(£>···<£>a<£>···®l with a acting on the i 's 
position of the tensor product. Notice that exp(iF n (a)) = 
exp(ia/V") 0n andy/n[F n {a),F n {b)]=F n ([a,b]). Then 

lim \\L(F n (a),F n (b))\\ =0. 

n — »oo 

The convergence is uniform over ||a||, ||6|| < C for some 
constant C . 

This Theorem is a key ingredient of the quantum cen- 
tral limit Theorem 36] and it is not surprising that it 
plays an important rolc in our quantum local asymp- 
totic normality result which is an extension of the lat- 
tcr. We apply the Theorem to two unitàries of the form 
U(u/x/n)® n = exp(i(u x a x + u v a y ))/^n)® n . We thus get 
information on the effect of the Uj(u/y/n) on the highest 
weight vectors of an irreducible representation. 

Corollary 5.3 For any unitary U and state r let 
Ad[í7](r) := UtU* and consider the rotated states 



r(u,v,j,n) := Ad 
r(u + v,j, n) := Ad 



Ui 



u 

y/n 

U + V 

y/n 



U, 



mm 



mm 



Then the following uniform convergence holds 
lim sup sup ||r(u, v, j,n) - r(u + v, j,n)\\i = 0. 

n ^°°u,ve/ 2 3'6J»,e 

Proof: By applying Theorem 15.21 to U(u/^fn)® n and 
U(v/s/n)® n , the first term of L{F n {a)) is 

/ \ SJn / \ <gm 
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and the second is 



U 2 = U [ — — ) exp 1 



2YH 



with 



[u, v] := [u x a x +u y a y ,v x a x + v y (7 y ] = 2{u x v y — U y v x )a z . 



The norm one distance between these two operators is 
going to as n is going to infinity, uniformly on (u, v) € 
I 2 . We may apply these operators on any pure state of 
(C 2 )®", in particular on for any j 6 J n ,t after 

block-diagonalization and preserve the uniform limit 



|Ad [^(lijXüD-Ad^diiXüDii 



0. (5.2) 



Now the action of exp í ^2^n^ J on 0' Ü ^ s s i m ply 
identity because \jj) is an eignevector of J JiZ . Thus 

Ad [U 2 ](\jj){jj\) =T(u + v,j,n). 

Togheter with (|5.2|) . this ends the proof. 

□ 

The following Lemma is a slight strengthening of a 
theorem by Hayashi and Matsumoto p9j . 

Lemma 5.4 The uniform convergence holds 



lim sup sup 



-|\/2/í - la u ) 



0. 



where |z) denotes a coherent state of the oscillator, and 
Q' u := (— 1% + iu x ) . Moreover for any sequence j n — ► oo 
we have 



lim K^V^-^II^ 



o. 



(5.3) 



The convergence holds uniformly over all sequences j n 
such that j n /n > c for some fixed constant c > 0, so in 
particular for j n S J n .t- 

Proof. We first prové the easier relation (|5.3(l . As both 
density matrices are diagonal we get 



\y. p o v* - ó°\\ - ^"^V / 

" í n 



2j„ 



fe=0 



fc=2j„+l 



as n — > oo. 

As for the first relation, let us denote |u,j, n) 
VjUj(-^)\jJ), then by and (fÏÏTj) we have 



;/-■ u..y. ») = w ' ( ^ )(sin(|u|A^)e^) fc (cos(|u|A/ïï)) 2 -' ' 



Now, the following asymptotical relations hold uniformly 
over j S J n ,e ■ 



sm 



u 



n / 

2j-k 



Jn 



(l + Oflufn" 1 )), 



exp(-fciK)(l + 0(|u|^)), 



((2/x-l)n) 
fe! 



-(l + 0(n" e )), 



and thus the coefficients connverge uniformly to those of 
the corresponding coherent states as n — > oo 



exp 



(2 M -l)|u| 2 \ (e^lulv^T) 



□ 

Proof of Provosition ïh^Ti The main idea is to notice that 
0° is a thermal equilibrium state of the oscillator and 
can be generated as a mixture of coherent states with 
centered Gaussian distribution over the displacements: 



1 



V2 



7TS 



\z)(z\d 2 z. 



(5.4) 



The easiest way to see this is to think of the oscillator 
states in terms of their Wigner functions. Indeed, the 
Wigner function of a coherent state is 



W z {q,p) = exp (-(q - \/2Rcz) 2 - (p - \/21mz) 



and thus the state given by (|5.4|) has Wigner function 
which is the convolution of two centered Gaussians which 
is again a centered Gaussian with variance equal to the 
sum of their variances 2s 2 + 1/2 which is equal to the 
variance of (jP for s 2 :=p/(2(l —p)). Similarly, 



1 



-|z-V2/í-1q u | 2 /2s 



^ 27TS 2 

Let us first remark that 
II^I^-^llx< 



(|z)(z|)d 2 z. (5.5) 



\pt-vm\\ 1 + 



Pi. 



i p, 



in Ml ) 



where Pj n — Vj n VT is the projection onto the image of 



Vj n , and 



lim sup sup W~Pj (/) U P 3 |L =0, 



because j n — > oo uniformly and Pj n converges to the iden- 
tity in strong operator topology (a tightness property). 
Thus it is enough to show that 

lim sup sup Lo!? — V* 4> u Vj II =0. 
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Now 



Ad 






M 








K 






Ad 


u jn ( 











< 



The first term on the right side of the inequality converges 
to zero by Lemma |5.4I uniformly for any sequence (j n ) 
such that j n £ J n ,í and does not depend on u. Using 
(|5.4|) and (|5.5|) we bound the second term by 



i r e -i Z i 2 /^|| A(U;ZiJn) n ld 2 z 

sy/zn J 

where the operator A(u, z,j n ) is given by 



A(u,z,j„) := Ad 



U, 



u 

i/n 



(yi\*){z\v jn ) 



y/2p, - ía u ) ( z + y/2p- la u 



V,, 



We analyze the expression under the integral. Let z G 
be such that = z/y/2p, — 1, then 



Ad 
Ad 



V jn Ad 



V jn Ad 



(V*Jz){z\Vj n ) - V*Jz + y/2p~ïa u )(z + y/2p~J ail \V j7 



< 



Ui 



u. 



u, 



z 

y/n 
u + z 

y/n 



Qjnjn) 0'nÍnl) ~ Ad 
(|jnJn)OnJn|)^:-|z)(z| 



(|jnjn>0'njn|) 



(lÍnÍn)O'njnl) V£ - |z + - la u )(z + ^2/i - la u | 



By Corollary 15.31 the first term on the right side con- 
verges to zero uniformly in (u, j n ) £ I 2 x J U:( . By Lemma 
I5.4l we have that the last two terms converge to zero uni- 
formly in (u, j n ) £ I 2 x J n f . Thus if we denote 

F n (z) := sup sup ||A(u,z, j n )\\ 1 

then < F n (z) < 2, lim, woo F„(z) = for all z £ R 2 , 
and by the Lebesgue dominated convergence theorem we 
get 

lim —L· / e -l-l 2 / 2s V„(z)d 2 z = 0. 

n->oo S ^27T J 

This implies the statement of the Proposition l5.ll 

□ 



6. CONSTRUCTION OF THE INVERSE 
CHANNEL S n 

To complete our proof of asymptotic equivalence as 
defined by l|1.4l) , we must now exhibit the inverse channcl 
S n which maps the displaced thermal states </) u of the 
oscillator into approximations of the rotated spin states. 
As the latter are block diagonal with weights p n (j) as 
defined in equation (|4.7|l , it is natural to look for S n of 



the form 

n/2 

Sn(4>)= PnV)S' n (tt>)®—, 

where S n are channels with outputs in Tíj. Moreover 
because Vj : Tíj — > L 2 (R) is an isometry we can choose 
S n such that 

SÍ(V jP V*)=p, (6.1) 

for all density matrices p on Tíj. This property does 
not fix the channel completely but it is sufficient for our 
purposes. Basically what we want is an inverse of the 
embedding Vj ■ V* used for the direct channel and one 
way to get this is as follows First block diagonalize 4> to 
get Pj4>Pj + Pj(f>Pj where Pj is the projection onto the 
image of Vj, i.e. Pj = VjV* , and note that this is a trace 
preserving completely positive map. This block diagonal 
state can be now seen as a state on the direct sum àlgebra 
6(P JJ L 2 (K))©6(P j L L 2 (]R)) and can be mapped to a state 

on B(Ti 3 ) by the channel V* ■ Vj © Sj- with Sj- arbitrary 
on the 'upper block'. The resulting composition is the 
channel S n satisfying property i|(j.lf) . 

Theorem 6.1 The following holds 

lim sup ||S„(> u )-p"l|i = 0. 
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Proof. As both p" and S n (cj> u ) are block-diagonal we may 
decompose their distance as 

n/2 

\\S n (r) - Pnh = E P«Ü)II^(0 U )-P"„I|1 
+ E PnÜ')ll^(^P°,nl7)-P°,nlll 

<2 e pn(i)+ e pnimr-VjPi^iu, 

where we have used at the last line that 5*^ is a contrac- 
tion and property (|6.1|l of S J n . Now the first sum is going 
to by H4.10fl and the second sum is also uniformly going 
to by use of Proposition l5.ll 

□ 



7. APPLICATIONS 

A. The optimal Bayes measurement is also 
asymptotically local minimax 

In this subsection we will introduce some ideas from 
the pointwise approach to state estimation. We show 
that the measurement which is known to be optimal for 
a uniform prior in the Bayesian set-up, is also asymptot- 
ically optimal in the pointwise sense. 

Using the jargon of mathematical statistics, we will 
call quantum statistical experiment (model) 0] a family 
{p e 6 M(C d ) : 9 e 9} of dcnsity matrices indexed by 
a parameter belonging to a set O. The main examples 
of quantum statistical experiments considered so far are 
that of n idcntical qubits 



A natural choice of distance between density matrices is 
related to the fidelity squarc 



F(p,af= [Tr ((^Vp) 1 



/2 



T := {p® n :pe AÍ(C 2 )}, 



the local model 

and its "limit" model 

Q 1 :={4> u ,uel 2 }, 

where / = [— a, a] , and p" and </> u are defined by (|l.l(l and 
(|1.2|1 . More generally we can replace the square I 2 by an 
arbitrary region K in the parameter space and obtain: 

Ç K :={f,ueíícl 2 }. 

We shall also make use of 

5:={f,uel 2 }. 



which is locally quadratic in first order approximation, 
i.e. 

i-nc^) 2 «-iiu-v|| 2 . 

n 

As we expect that reasonable estimators are in a local 
neighborhood of the true state we will replace the fidelity 
square by the local distance 

d(u, ü) = ||ü — u|| 2 . 

and define the risk of a measurement-estimator pair as 
Rm(u, ü) = (d(u, ü)), keeping in mind the factor 1/n re- 
lating the risks expressed in local and global parameters. 

Similarly to the Bayesian approach, we are interested 
in estimators which have small risk everywhere in the 
parameter space and we define a worst case figure of mèrit 
called minimax risk. 

Definition 7.1 The minimax risk of a quantum statis- 
tical experiment £ over the parameter space for loss 
function d{9 1 0), is defined as 



C{£) = inf supR M (0,9). 

M,è eee 



(7.1) 



where the infïmum is taken over all measurement- 
estimator pairs {M,9), and Rm(9,9) — {d{9,9)). 

The minimax risk telis us how difficult is the model and 
thus we expect that if two models are "close" to each 
other then their minimax risks are almost equal. The 
"statistical distance" between quantum experiments is 
defined in a natural way with direct physical interpre- 
tation and such a problem has been already addressed 
in |43l | for the case of a quantum statistical experiment 
consisting of a finite family of pure states. 

Definition 7.2 Let £ = {p e e M(C d ) : 9 € 0} and 
T = {t 6 E M(C p ) : 9 6 0} be two quantum statistical 
experiments (models) with the same parameter space 0. 
We define the deficiències 

ó(£,F)=mïsup\\T(p e )-T e \\ 1 , 
1 eee 

5{T,£) =infsu P ||//-S(T 9 )|| 1 , 
s eee 

where the infimum is taken over all trace preserving chan- 
nels T : M(C d ) -> M(C P ) and S : M(C P ) -> M(C d ). 

With this terminology, our main result states that for 
any bounded interval /: 

lim max^,^),^,^)) = 0. (7.2) 
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As suggested above, the deficiency has a direct statisti- 
cal interpretation: if we want to estimate 9 in both sta- 
tistical experiments £ and T and we choose a bounded 
loss function d(9, 9) < K then for any measurement and 
estimator 9 for T with risk Rm(9,9) = (d(9,9)) we 
can find a measurement N on £ whose risk is at most 
Rm(9,9) + KS(£,!F). Indeed if we choose T such that 
the infimum in the definition of 8(£, T) is achieved, we 
can map the state p s through the channel T and then 
perform M to obtain an estimator 9 for which 



R N (9, 9) = (d(9, 9)) = J d{9, 9)Tt ( T(p")M(d9) ) < 



d(9,9)Tr t M(d9) + NU|ïV) - r w ||i < 



R M (9,d) + K6(£,F). 

This means that the difhculty of estimating the param- 
eter 9 in the two models is comparable within a factor 
S(£,J-). With the above definition of the minimax risk 
and using the convergence 1)7.2(1 we obtain the following 
lemma. 

Lemma 7.3 Let I = [—a, a] with < a < oo, then 
hm C(F I n ) = C(G I ). 

n — >oo 

The minimax risk for the local family is a figure of 
mèrit for the "local difficulty" of the global model T n . It 
converges asymptotically to the minimax risk of a family 
of thermal states. However this quantity depends on the 
arbitrary parameter / = [— a, a] which we would like to 
remove as our last step in defining the local asymptotic 
minimax risk 

a.a.m.(^n : n G N) := lim lim C(J^) = lim C(Q T ). 

a — s-oo n — >oo a — >oo 

This quantity depends in principle on the state which 
is at the center of the local neighborhood. However by 
invariance under rotations, the risk is constant for a given 
pair of eigenvalues of the density matrix. Now, as one 
might expect the minimax risks for the restricted families 
of thermal states will converge to that of the experiment 
with no restrictions on the paramaters. The proof of this 
fact is however non-trivial. 

Lemma 7.4 Let I = [~a,a], then we have 
lim CiG 1 ) = C(G) 

a — >oo 

Moreover the heterodyne measurement saturates C(G), 
and thus C(G) is equal to the Holevo bound. 

Proof. The inequality in one direction is easy. For any es- 
timator, sup uS/ 2 í?aí(u, ü) < sup u6R 2 í?m(u, ü), so that 
C-iG 1 ) < C{G) and the same holds for the limit. By 
the same reasoning, for any K\ C K% C R 2 we have 
C(G Kl ) < C(G K2 ). 



When calculating minimax bounds we are interested 
in the worst risk of estimators within some parameter 
region K, and this worst risk is obviously higher than 
the Bayes risk with respect to the probability distribution 
with constant density on K . We shall work on B(0, c + b) 
the ball of center and radius (c + b), with b > c, and 
denote our measurement M with density m(ü)dü. In 
general M need not have a density, but this will ease 
notations. Then 



sup í?m(u, ü) > 

uGB(0,c+6) 

du dü 



u-ü|| 2 Tr(çi u m(ü)). (7.3) 



'B(0,c+b)xK 2 n i c + 

We fix the following notations 

/(£>)= / dudv||x-y|| 2 Tr(,£ u m(v)), 
Jv 

g(V) = / dudvTr (çi u m(v)) , 
Jv 

and define the domains 

Oi = {(u, ü)|u e B(0, c + b),üe R 2 } 
V 2 = {(u + k, k)|u G B(0, c), k G B(0, b)} 
2?3 = {(u,u + h)|u€fl(0,6-c),heB(0,c)} 
Vi = {(u, u + h)|u G B(0, b - c), h G R 2 \B(0, c)}. 



Notice the following relations: 

V3C.V2C.Vi, ViCVi\V 2 . 
Then (|7.3|) can be rewritten as 

1 



(7.4) 



SUp Í?aí(u, Ü) > 

uS B(0, c+b) 



n(b + c)'- 



The following inequalities follow directly from the def- 
initions: 

f(V 2 ) < c 2 g(V 2 ) f(V 3 ) < c 2 g(V 3 ) 

f(Pi) > c 2 g(V±) g(V 4 ) + g(V 3 ) = n(b - cf . 

Using these and l|7.4|) . we may write: 

1 s 1 



7r(c + Ò) 2 
1 

> 



f(Vi) > 



ir(c + b) : 



(f(V 2 ) + f(V A )) 



(f(V 2 ) + c 2 g(V 4 )) 



7r(c + Ò) 2 

(ò-c) 2 f f(V 2 ) g(V 2 ) 2 _ 2 g(V 3 ) 
(b + c) 2 \g(V 2 ) ir(b - c) 2 C tt(ò-c) 2 



> 



> 



(b~c) 2 , , 
(b + cf 

(b-c) 2 f(V 2 ) 
(b + cYg(V 2 )' 



g(V 3 ) ff(V 2 ) 



tt(6-c) 2 \g(V 2 ) 



(7.5) 
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We analyze now the expression f (V2) / g{T>2) ■ By using 
the definition (|1.2I) of the displaced thermal states (j> u we 



get that Tr 



AU + k, 



t(l)] = Tr |> k m u (i)]. where 



m u (l) := D{-^2ii - la ví )m{\)D{^2 l· L - ía u ). 



Then 



dudkTÏ 



I B(0,c)xB(Q,b) 

where we have written 



i(k)] = Tr (p c m b 



Mu, 



B(0,c) 



Upon writing v c := J 1 



rrii, = / mk(k)dk. 

'B(0,b) 

u|| 2 u du, we get similar ly 



B(0,c) 

/(P2) — Tr [w c mf,]. Note that by rotational symmetry v c 
and 4> c are diagonal in the number operator eigenbasis, 
so without restricting the generality we may assume that 
ràfc is also diagonal in that basis: rhb — ^2^Pk\k){k\. We 
have then 

/(Pa) = E k ^Pk(k\vc\k) > . nf (fc|^|fc) 
E* 6 NP*M&l fc > -*6N(fe|0 c |fe)- 

The infimum on the right side is achieved by the vac- 
uum vector. By Lemma [7.51 this fact follows from the 
inequality 



(k\cj) Ul \k) (O|0 Ul |O) 

(k\<t> n *\k) ~ (o\4> u *\o)' 



l U l|| > ll U 2| 



which can be checked by explícit calculations. 

Letting now c and b go to infinity with c = o(b) and 
using l|7.5|) . we obtain that 



lim C(Ç a ) > 



/ R2 (oi0"lo)llu|| 2 du 
; R2 (o|^|o)du ' 



which is exactly the pointwisc risk of the heterodyne mea- 
surement H(du) = h(u)du whose density is 

h(u) = (2/i- l)^(-y2^Ta u )|0)(0|^(-v^7^Ta u ). 

By symmetry this pointwise risk does not depend on the 
point, so that C{G) < Rh(vl, ü). And we have our second 
inequality: lim^oc C(G a ) > C{Q). 

Moreover, the heterodyne measurement is known to 
saturate the Holevo bound for G — Id and the Cramér- 
Rao bound for locally unbiased estimators |29l l39| . We 
conclude that the local minimax risk for qubits is equal to 
the minimax risk for the limit Gaussian quantum experi- 
ment which is achieved by the heterodyne measurement. 

□ 

Lemma 7.5 Let p and q be two probability densities on 
[0, 1] and assume that 



p(xi) > q(xi) 
p(x 2 ) ~ q{x 2 y 

Then J x 2 p{x)dx > J x 2 q(x)dx. 



x\ > x 2 - 



Proof. It is enough to show that there exists a point 
xq 6 [0, 1] such that p{x) < q{x) for x < xq and p{x) > 
q(x) for x > xç). Now, if p{x) < q(x) then by using 
the assumption we get that p(y) < q(y) for all y < x. 
Similarly, if p{x) > q{x) then p{y) < q{y) for all y > x. 
This implies the existence of the crossing point Xq. 

□ 

By putting the last two lemmas together we obtain the 
following. 

Proposition 7.6 The local asymptotic minimax risk 
Ci.a.m^Fn ■ fi € N) for the qubit state estimation problem 
is equal to the minimax risk C(G) for the corresponding 
quantum Gaussian shift experiment which is achieved by 
the heterodyne measurement. 

The natural question is now the following: is there a 
sequence of measurement-estimator pairs for the qubits 
which achieves this the risk Ci. a . m (J r n : n € N) asymp- 
totically for all local neighborhoods simultaneously, i.e. 
without prior knowledge of the center p° of the í/y/n ball 
within which the true state lies. Intuitively, the following 
procedure seems natural: use the local asymptotic nor- 
mality to transfer the heterodyne measurement from the 
space of the oscillator to that of the qubits and in this 
way achieve the desired asymptotic risk. However this re- 
quires the knowledge of the local neighborhood on which 
the convergence holds. In order to obtain this informa- 
tion about the state we need to 'localize' the state by per- 
forming a first stage of (rough) measurements on a small 
proportion of the systems of order o(n) and then per- 
forin the (optimal) heterodyne type measurements cor- 
responding to the local neighborhood of the first stage 
estimator. In order to make this argument rigurous we 
need some finer estimates on the region in which local 
asymptotic normality holds and we leave this problcm 
for a separate work. 

However, there exists another measurement which 
achieves the risk Cí. a . m (J-" n : n € N), namely the opti- 
mal measurement from the Bayesian point of view dis- 
cussed in pü. |2^ | . The connection between the local and 
Bayesian approaches is discussed in more details in the 
next Subsection to which we refer for the appropriatc 
dcfmitions. In particular, the next proposition can be 
better understood after reading the next Subsection but 
we state it here because it is a direct consequence of the 
results derived in this Subsection. 

Let us denote by (M n ,p n ) the measurement-estimator 
pair from 0) HH which are optimal from the Bayesian 
point of view. 

Proposition 7.7 The optimal measurement-estimator 
pair (M„,ü n ) in the Bayesian setup is an asymptotically 
local minimax estimation scheme. That is for any p 

lim nR Mn {PiPn) = C'i.a.m^n ! 71 G N), 

n — >oo 

where R,M n (p, p n ) is the risk with respect to teh fldelity 
distance. 
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Proof. The pointwise risk of (M n ,p n ) is known to con- 
verge to that of the heterodyne measurement [3 ■ The 
rest follows from Lemma 17.31 and Lemma IT4l 

□ 

B. Local asymptotic equivalence of the optimal 
Bayesian measurement and the heterodyne 
measurement 

In this subsection we will continue our comparison of 
the pointwise (local) point of view with the global one 
used in the Bayesian approach. The result is that the 
optimal SE/ (2) covariant measurement j^iH^ converges 
locally to the optimal measurement for the family of dis- 
placed Gaussian states which is a heterodyne measure- 
ment 39] . Together with the resulta on the asymptotic 
local minimax optimality of this measurement, this closes 
a circle of ideas relating the different optimality notions 
and the relations between the optimal measurements. 

Let us recali what are the ingredients of the state es- 
timation problem in the Bayesian framework |24j . We 
choose as cost function the fidelity squared F(p, a) 2 = 
Tr(^/ y/pcTyfp) 2 and fix a prior prior distribution tt over 
all states in C 2 which is invariant under the SE/ (2) sym- 
metry group. Givcn n idcntical systems p®" we would 
like to hnd a measurement M n - whose outcome is the 
estimator p n - which maximizes 

Rn,n ■= j {F( Pn ,p) 2 )n{dp). 

By the SU (2) invariance of tt, the optimal measurement 
can be chosen to be SU (2) covariant i.c. 

UM n (da)U* = M n (U*daU), 

and can be described as follows. First we use the decom- 
position (|4.5(l to make a "which block" measurement and 
obtain a result j and the conditional state p^ n as in l|4.ti[) . 
This part will provide us the eigenvalues of the estimator. 
Next we perform block-wise the covariant measurement 
Mj,n(d~s) = mj_ n (~s )d~s with 

m^Çt) := (2j + vujrfyiïïWjCf) ® 1, 

whose result is a unit vector "s* where U(~s) is a uni- 
tary rotating the vector state \~s) to | f"). The complete 
estimator is then p n = |(1 + ^~s~a). 

We pass now to the description of the heterodyne mea- 
surement for the quantum harmònic oscillator. This mea- 
surement has outeomes u 6 M 2 and is covariant with 
respect to the translations induced by the displacement 
operators D(z) such that H(du) = h(u)du with 

h(u) := (2p - í)D{~^2p~ía u )\0)(0\D(^2p~ía u ). 

Using Theorem 11.11 we can map H into a measurement 
on the n-spin system as follows: first we perform the 



which block step as in the case of the SE/(2)-covariant 
measurements. Then we map p^ n into an oscillator state 
using the isometry Vj (see l|5.1|l ). and subsequently we 
perform H. The result u will dehne our estimator for 
the local state, i.e. 

We denote by H n the resulting measurement with vàlues 
in the states on C 2 . 

The next Theorem shows that in a local neighborhood 
of a fixed state p°, the SE/(2)-covariant measurement M n 
and the heterodyne type measurement H n are asymptot- 
ically equivalent in the sense that the probability dis- 
tributions P(M n ,p) and P(H n , p) are close to each other 
uniformly over all local states p such that \\p — p°\\ i < 
for a fixed but arbitrary constant C < oo. 

Theorem 7.8 Let p° be as in 1)4. l[l . and let 

B n (l) = {p v ^ l :vel 2 }, , | J| < oo 

be a local family of states around p° . Then 

lim sup \\P(M n ,p)-P(H n ,p)\\ 1 = 

n ^°° p£B n (I) 

Proof. Note first that both P(M n ,p) and P(H n ,p) are 
distributions over the Bloch sphere and the marginals 
over the length of the Bloch vectors are identical because 
by construction the first step of both measurements is 
the same. Then 

\\P{M n ,p)-P{H n ,p)\\ x = 

^PnÜ) / \Tr(p jin (m jin (~s) - h jtn Çs)))\d~s. 

3 J 

According to H4.1U|) we can restrict the summation to 
the interval J n ^ around j = n(p — |). By Theorem 
11.11 we can replace (whenever needed) the local states 

by their limits in the oscillator space </> v with an 
asymptotically vanishing error, uniformly over v G I 2 . 

We make now the change of variable s — > u where u € 
R 2 belongs to the ball |u| < 2-y/wr, and is the smallest 

vector such that U — UÇ~s). 

The density of the SU (2) estimator with respect to the 
measure du is 

where J is the determinant of a Jacobian related with 
the change of variables such that J(0) = 1. 

Similarly the density of the homodyne-type estimator 
becomes 

h jtn {u) := V h ( u + 2k ^^) V i I J fc.»( u )l> 
keíi ^ ' 
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because displacements in the same direction which dif- 
fer by múltiples of 2y/n-K lead to the same unitary on 
the qubits. Here Jfc jrl (u) is again the determinant of the 
Jacobian of the map from the fc-th ring to the disk, in 
particular Jo, n (u) = 1. 

The integral becomes then 



|u|<27rVn 



Tr 



{Pj!rf í { m i,n{ u ) - hj, n (u)) 



du. 



We bound this integral by the sum of two terms, the 
first one being 



|u|<27r v /rï 



Tr 



/ v/y'n/ 



(u)-^(u)) 



du, 



where hj(u) is just the term with k = in hj tn (u). By 
Lemma I5.4I for any fixed u we have mj in (u) — > h(u) 
uniformly over j £ 3 n .<L- Using similar estimates as in 
Lemma |5 .4! it can be shown that the function under the 
integral is bounded by a fixed integrable function <?(u) 
uniformly over ve/ 2 , and then we can use dominated 
convergence to conclude that the integral converges to 
uniformly over v € I 2 and j £ J n ,e- 
The second integral is 



/ 



'|u|<2·7r v /rï 

which is smaller than 



|u|>27r v /rï 



Tr( P ;//%(u)-%„(u)) 



du, 



Tr 



(P%^h(u) 



du, 



which converges uniformly to 0. This can be seen by 
replacing the states with </> v which are "confined" to a 
fixed region of the size I 2 in the phase space, while the 
terms h(u) are Gaussians located at distance at least 
2i:^/n from the origin. 

Putting these two estimates together we obtain the 
desired result. 

□ 

Remark. The result in the above theorem holds more 
generally for all states in a local neighborhood of p° but 
for the proof we need a slightly more general version of 
Theorem II . II where the eigenvalues of the density matri- 
ces are not fixed but allowed to vary in a local neigh- 
borhood of (/i, 1 — p,). This result will be presented in a 
future work concerning the general case of d-dimensional 
states. 



C. Discrimination of states 

Another illustration of the local asymptotic normality 
Theorem is the problem of discriminating between two 
states p + and p~ . When the two states are fixed, this 
problem has been solved by Helstrom 0| , and if we are 
given n systems in state p®" then the probability of error 



converge to exponentially. Here we consider the prob- 
lem of distinguishing between two states p^ which ap- 
proach each other as n — > oo with rate ||p£ — Pü\\i ~ 
In this case the probability of error does not go to 
because the problem becomes more difficult as we have 
more systems, and converges to the limit problem of dis- 
tinguishing between two fixed Gaussian states of a quan- 
tum oscillator. 

This problem is interesting for several reasons. Firstly 
it shows that the convergence in Theorem ll . ll can be used 
for finding asymptotically optimal procedures for various 
statistical problems such as that of parameter estimation 
and hypothesis testing. Secondly, for any fixed n the opti- 
mal discrimination is performed by a rather complicated 
joint measurement and the hope is that the asymptotic 
problem of discriminating between two Gaussian states 
may provide a more realistic measurement which can be 
implcmented in the lab. Thirdly, this example shows 
that a non-commuting onc-parameter familics of states 
is not "clàssica!' as it is sometimes argued, but should 
be considered as a quantum "resource" which cannot be 
transformed into a classical one without loss of informa- 
tion. More explicitly, the optimal measurement for esti- 
mating the parameter is not optimal for other statistical 
problems such as the one considered here, and thus dif- 
ferent statistical decision problems are accompanied by 
mutually incompatible optimal measurements. 

Let is recali the framework of quantum hypothesis test- 
ing for two states pr : we consider two-outcomes POVM's 
M = (M-,M+) with < M + < 1 and M- = 1 - M+ 
such that the probability of error when the state is p~ is 
given by Tr(Aí + p~),and similarly for p + . As we do not 
know the state, we want to minimize our worst-case prob- 
ability error. Our figure of mèrit (the lower, the better) 
is therefore: 



R(p + ,p~) = infmax{Tr(p + M_),Tr(p+M_)} 



Now we are interested in the case when p ± — p„ as 
defined in (|l.lf) . and in the limit p± — cf> ±u (recali that 
both and 4> u depend on p). We then have: 

Theorem 7.9 The following limit holds 

hm R(p»,p-») = R(^A- U )- 

n — >oo 

Moreover for pure states this limit is equal to 
( 1 — (1 — e -4 !"! ) 1 / 2 ^ /2 which is strictly smaller than 

1/2 — er/(|u|) which is the limit if we do not use collective 
measurements on the qubits. Here we have used this con- 
vention for the error function: erf(x) = e~ l /y/ndt. 

Proof. Let M be the optimal discrimination procedure 
<f> ±u . Then we use the channel T n to send p„ to 
states of the oscillator and then perform the measure- 
ment M. By Theorem OI ||0 ±u - T„(p^ u )||i -> so 



that Tr(T„(p± u )M T ) -> Tr ( 
asymptotically optimal for 



M-, 



Thus M o T„ is 



±u 
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Now for pure states and |0_) the optimal mea- 
surement is well-known |45l l46l| . It is unique on the span 
of these pure states and arbitrary on the orthogonal. If 
we choose the phase such that (ip-\i/}+) > 0, then M + is 
the projector on the vector 

10+) + \P-) + 10+) - kM 

2 x /l+(V>_|</>+) 2^1-^-1^+) 
and the associated risk is 

Vl-|(^+|V^)| 2 ) 

Now in our case, in the limit experiment, U is the co- 
herent state |^ u ) = e-H 2 / 2 \u\ n /Vriï \n). So that 

z — ' 7l! 

n 

and i?(0 u , 0- u ) = è (l - Vi -e~ 4 M^ . 

We would like to insist here that the best measurement 
for discrimination is not measuring the positive part of 
the position observable Q (we assume by symmetry that 
±u is on the first coordinate), as one might expect from 
the analogy with the classical problem. Indeed if we 
meausure Q then we obtain a classical Gaussian variable 
with density p(x) — e - ^ - ") j ^pn and the best guess at 
the sign ± has in this case the risk 1/2 — er/(|u|). 

This may be a bit surprising considering that mea- 
suring Q preserves the quantum Fisher information. The 
conclusion is simply that the quantum Fisher information 
is not an exhaustive indicator of the statistical informa- 
tion in a family of states, as it may remain unchanged 
even when there is a clear degradation in the inference 
power. This example fits in a more general framework of 
a theory of quantum statistical experiments and quantum 
decisions |47| . 

□ 



states can be prepared deterministically by using contin- 
uous time measurements performed in the environmcnt 
and real time feedback on the spins. Without going in 
the details, the bàsic idea is to describe the evolution of 
identically prepared spins by passing first to the coherent 
state picture. There one can easily solve the stochastic 
Schròdinger equation describing the evolution (quantum 
trajectory) of the quantum oscillator conditioned on the 
continuous signal of the measurement device. The solu- 
tion is a Gaussian state whose center evolves stochasti- 
cally while one of the quadratures gets more and more 
squeezed as one obtains more information through the 
measurement. Using feedback one can then stabilize the 
center of the state around a fixed point. 

This description is of course approximative and holds 
as long as the errors in identifying the spins with Gaus- 
sian states are not significant. The framework developed 
in the proof of Theorem 11.11 can then be used to make 
more precise statements about the validity of the rcsults, 
including the squeezing process. 

Perhaps more interesting for quantum estimation, such 
measurements may be used to perform optimal estima- 
tion of spin states. The idea would be to first localize the 
state in a small region by performing a weak measure- 
ment and then in a second stage one performs a hetero- 
dyne type measurement after rotating the spins so that 
they point approximately in the z direction. We believe 
that this type of procedure has better chances of being 
implemented in practice compared with the abstract co- 
variant measurement of |24ll29j|. 
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D. Spin squeezed states and continuous time 
measurements 



In an emblemàtic experiment for the field of quantum 
filtering and control [35( it is shown how spin squeezed 
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